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We propose a new theory of massive gravity with only two propagating degrees of freedom. While 
the homogeneous and isotropic background cosmology and the tensor linear perturbations around 
it are described by exactly the same equations as those in the de Rham-Gabadadze-Tolley (dRGT) 
massive gravity, the scalar and vector gravitational degrees of freedom are absent in the new theory 
at the fully nonlinear level. Hence the new theory provides a stable nonlinear completion of the self- 
accelerating cosmological solution that was originally found in the dRGT theory. The cosmological 
solution in the other branch, often called the normal branch, is also rendered stable in the new 
theory and, for the first time, makes it possible to realize an effective equation-of-state parameter 
different from (either larger or smaller than) —1 without introducing any extra degrees of freedom. 


Introduction. — Since the seminal work by Fierz and 
Pauli [1], especially in the recent years, much theoreti¬ 
cal effort in cosmology has been put in order to develop 
theories of massive gravity [2, 3|. These theories were 
indeed able to introduce, at non-linear level, the desired 
five modes necessary to describe a massive graviton in a 
Lorentz invariant way. In other words, these theories are 
free from the so called Boulware-Deser ghost [4|, which 
had been thought to plague any theories of massive grav¬ 
ity. Together with this first success, much work came in 
order to see whether these same theories could be viable. 
Unfortunately, these theories suffer from instability on 
some key backgrounds, such as the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) universe [j|. In this regard 
several attempts have been analyzed to find stable cos¬ 
mological solutions in massive gravity: 1) abandoning the 
homogeneity and/or the isotropy of cosmological models; 
2) changing the theory by introducing new fields inter¬ 
acting with gravity in a way as to save the theory. It 
proved difficult, even with these attempts, to find a the¬ 
ory of massive gravity with a theoretically consistent and 
experimentally viable cosmology. 

In this letter, we present a new theory of massive grav¬ 
ity which modifies general relativity in a minimal way. 
We will perform this by looking for a theory with only two 
tensor modes, which are massive. This will make FLRW 
backgrounds (including de Sitter) stable and viable as in 
the standard cosmology. Indeed the tensor modes of the 
gravity sector will be massive, whereas there are no scalar 
and vector propagating modes in the gravity sector. In 
order to achieve this goal we will not impose the Lorentz 
symmetry, so that a massive graviton does not need to 
have five degrees of freedom any longer. 

Precursor Theory. — In order to define the theory 
we will make use of the lapse N, shift N‘\ and the three- 
dimensional vielbein e 1 j as basic variables. We can then 
introduce the three-dimensional metric as 

Tij — &IJ & i & j • (1) 

Hereafter, I,j£ (1, 2, 3} so as i and j. Out of the vari¬ 


ables introduced so far, we can build a four-dimensional 
vielbein as 


le A II - 


N 0 T 

eW e J i 


and a four-dimensional metric as 
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where '!) A b is the Minkowski metric tensor, so that 

S oo = -^ 2 +7#W, 

goi = UijN 3 = g i0 , gij = t ij , (4) 

corresponding to the metric tensor written in the ADM 
variables. We also introduce a non-dynamical four¬ 
dimensional vielbein built out of a non-dynamical lapse 
M, a non-dynamical shift M l , and a non-dynamical 
three-dimensional vielbein E 1 j , as follows: 


( M 0 T \ 
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The four-dimensional vielbein of the form often 
called the ADM vielbein, has 13 independent compo¬ 
nents, as opposed to 16 independent components of a 
completely general vielbein in four-dimensions. The 
missing 3 components are the boost parameters that 
would transform the vielbein of the form m to a gen¬ 
eral vielbein. Therefore, by choosing the form (O for the 
vielbein, we introduce a preferred frame and thus explic¬ 
itly break the local Lorentz symmetry. 

We now introduce a precursor action, which will then 
be used as the starting point to define the theory: 


Ve = ~U~ f d i x y /~^gIZ[g l 
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+ je ABCV e a ^ s E A a E B 0 e c 1 e D s 
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where H[g^] is the four-dimensional Ricci scalar for the 
metric g^ v and the Levi-Civita symbol is normalized as 
eoi 23 = 1 = —e 0123 . The precursor action would be ex¬ 
actly the same as that for the dRGT massive gravity if 
e^a were a general, i.e. totally unconstrained, vielbein in 
four-dimensions. At the level of the definition of the pre¬ 
cursor theory, the only difference from the dRGT theory 
is thus that the four-dimensional vielbein is restricted to 
the form (J 2 J). 

Having given the action for the precursor theory, it is 
straightforward to write down its Hamiltonian. The pre¬ 
cursor Hamiltonian turns out to be linear in N and N 1 
and independent of their time derivatives. One can thus 
safely consider N and N l as Lagrange multipliers, and 
the phase space to be considered here then consists of 
9 x 2 = 18 variables, e L k and their conjugate momenta 
H k l • The coefficients of N and N l define primary con¬ 
straints, that we denote as — Hq and —Hi, respectively. 
The rank of the 4x4 matrix made of Poisson brackets 
among them is two, leading to two secondary constraints, 
which we denote as C T (r = 1,2). Combined with other 
six (three primary and three secondary) constraints, that 
we name as , pl MN 1 and Y 1 MJV 1 , associated with a sym- 
metry condition on the vielbein e L k , it is deduced that 
the physical phase space is 18 — 6 — 6 = 6 dimensional 
and that the number of physical degrees of freedom in 
the precursor theory is three at the fully nonlinear level. 

Minimal Theory.— While the precursor theory it¬ 
self is interesting, we further proceed to remove one more 
degree of freedom to define a theory with only two de¬ 
grees of freedom, that we call the minimal theory of mas¬ 
sive gravity. From now on, we will fix the units so that 
Mp = 2 . Also, we neglect the entirely non-dynamical 
part proportional to cq. 

The minimal theory is defined in the Hamiltonian lan¬ 
guage by imposing four constraints, which we denote as 
Cq and Cj and are defined in (0 below, on the precursor 
theory. Only two combinations among these four con¬ 
straints are new since the other two independent combi¬ 
nations are C T ~ 0 (r = 1,2), that already exist in the 
precursor theory. Hence the Hamiltonian of the minimal 
theory is 

H = J d 3 x[-NH Q - N'‘H, + m 2 MHi 

+XC 0 + XCi + a MN V lMN] + PmnY^} , (7) 

where N, N l , A, A 1 , a mat (antisymmetric) and Pmn (an¬ 
tisymmetric) are 14 Lagrange multipliers. This is a con¬ 
strained version of the precursor Hamiltonian, because 


we have added two additional constraints. As a conse¬ 
quence, on the constrained surface the Hamiltonian den¬ 
sity reduces only to H ss H\ = f d 3 xm 2 MT-Li- Each 
constraint has a specific meaning. The following terms 
are all derived from the precursor theory, 


Hq = Hq R — 771 2 "Ho , 
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out of which the precursor Hamiltonian is H pre = 
f d 3 x[-NH Q - WHi + m 2 M'Hi + A T C T + a M NV [MN] + 
PmnY\- mn }]. Here, r = 1,2, Vj is the spatial co¬ 
variant derivative compatible with 7 ,j, det 7 ^, 

Vl = VfetJ-, lij = S IJ E I i E J j, it* = 5 IJ W ie j k , Wr 

is the canonical momentum conjugate to e 1 j , and 


Y I J = E I k e J k , and X/ = e z k E J k , ( 8 ) 


satisfying Yi L X l j = Sj. 

Throughout the present letter, for simplicity we adopt 
the unitary gauge so that M, M l E I ; and E 1 j are only 
functions of the coordinates. This makes Hq and "Hi 
explicitly time-dependent. The remaining constraints, Cq 
and Ci , are then defined as 

BH 

C 0 = {H 0 ,H 1 } + -^-, Ci = {Hi, Hi} . (9) 

The two constraints C T ~ 0 (r = 1,2) in the precur¬ 
sor theory can be written as linear combinations of these 
four constraints. Therefore, only the remaining two com¬ 
binations are new. In other words, the minimal theory is 
defined by adding two additional constraints to the pre¬ 
cursor theory. The set of two new constraints removes 
one degree of freedom from the precursor theory. Since 
the precursor theory has three degrees of freedom, this 
means that the minimal theory has only two degrees of 
freedom. 

Rigorously speaking, what we have proved here is that 
there are enough number of constraints, meaning the in¬ 
equality, (number of d.o.f.) < 2, holds. One might in fact 
worry that the consistency of the additional two con¬ 
straints with time evolution might lead to further sec¬ 
ondary constraints, overconstraining the theory. This is 
not the case since, as we shall see explicitly in the next 
section, there are two (and only two) propagating degrees 
of freedom around cosmological backgrounds, meaning 
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another inequality, (number of d.o.f.) > 2, holds. By 
combining the two inequalities, we thus have the equal¬ 
ity, (number of d.o.f.) = 2, holds. It is also possible to 
prove the absence of further secondary constraints, and 
thus the presence of two physical degrees of freedom, in 
a more formal way by calculating the determinant of the 
14 x 14 matrix-operator made of Poisson brackets among 
all the fourteen constraints 

Having defined the minimal theory with only two de¬ 
grees of freedom by its Hamiltonian, it is straightfor¬ 
ward to calculate the corresponding action via a Legen¬ 
dre transformation. It should be noticed that since in the 
constraints, e.g. Co, the canonical momenta are present, 
on integrating out, e.g. the auxiliary field A, the resulting 
Lagrangian (or Hamiltonian) would acquire a structure, 
which would have a kinetic structure for 7 ^ essentially 
different from the Einstein-Hilbert one. While this mod¬ 
ification is essential for the absence of helicity -0 and -1 
gravitational modes, the modification to the kinetic term 
for the gravitational waves is suppressed by m 2 /Mp and 
thus is negligible. A nontrivial point to be noticed in 
this regard is that the relation between the canonical 
momenta H+ and the time derivative of the spatial viel- 
bein e 1 j is modified due to the additional terms in the 
Hamiltonian [@|. One can analyze the behavior of the 
theory, e.g. the cosmological evolution, by using either 
the Hamiltonian or equivalently the Lagrangian. 

Cosmology and phenomenology. — Let us con¬ 
sider a simple case to show the behavior of the theory, 
namely a general homogeneous and isotropic (FLRW) 
cosmological background with flat spatial metric, driven 
by the graviton mass term and matter fields minimally 
coupled to the metric g^. It is rather easy to see that 
on this background, Ci ~ 0 are trivial (because of homo¬ 
geneity of the background) and Co ~ 0 leads to 

(c 3 + 2c 2 A + c\X 2 ){X + NHX - MH) = 0 , (10) 


where A' = a/a is the ratio of the scale factors of the 
three-dimensional metrics 77 ,- and 7 ^ respectively, and H 
is the Hubble expansion rate (not the Hamiltonian), i.e. 
H = a/ ( aN ). This is exactly the same as the well-known 
constraint equation obtained from the Bianchi identity in 
the dRGT theory. Two branches of solutions thus exist, 
corresponding to the two factors of the left hand side 
of (HUD- The self-accelerating branch is defined by those 
values of the constant A' which satisfy 


X = X± = 


-c 2 ± \Jc\ - cic 3 

Cl 


( 11 ) 


while the normal branch is defined by setting X+NHX— 
MH = 0. In order to determine the Lagrange multiplier 
A, we demand that 


dXt) 

dt 


{Tl 0 , H} + — Co + 


d 3 y\{Ko,Co(y)} + ..., 


should vanish, where {7?-o,Co} ^ 0 and the neglected part 
vanishes because of the symmetry of the background. On 
requiring ~ 0 and imposing Co ~ 0 , we then find that 
A k 0 on the background. The remaining independent 
equation is, after re-inserting units, 

3 M^H 2 = IH-Mx ( C4 + 3c3 x + 3c 2 AT 2 + ciA 3 ) + p, 

(! 2 ) 

where p is the energy density of matter minimally coupled 
to the metric g^. This is exactly the same as the Fried¬ 
mann equation in the dRGT theory. There is no other 
independent equation for the background, essentially be¬ 
cause both the Bianchi identity and the constraint Co ~ 0 
lead to the same equation m- Needless to say, when we 
constructed the minimal theory in the previous section, 
we carefully chose Co so that this is the case. 

While the homogeneous and isotropic cosmological 
background solutions are exactly the same as those in 
the dRGT theory, perturbations around them behave 
completely differently. It is known that, in the standard 
dRGT theory, all homogeneous and isotropic background 
in both branches are plagued by ghosts, either in the 
helicity-0 or helicity-1 sector, and thus unstable [§]. On 
the contrary, in the new theory there is no physical degree 
of freedom in the helicity -0 or helicity -1 sector (namely 
no scalar or vector modes, according to the standard 1+3 
decomposition of the perturbation for the metric tensor) 
and thus those (would-be) ghosts are absent, rendering 
the cosmological background absolutely stable. In the 
minimal theory, the only existing two propagating modes 
reduce to the tensor modes, whose quadratic action can 
be written as 


S = 




d 4 xNa 3 


A=+,x' 


h{ (< dh x ) 2 


N 2 a 2 


where 


p 2 = 77 ra 2 A (c 2 A + c 3 ) + (ciA + c 2 ) 
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(13) 

(14) 


gives the mass p of the tensor modes on this background 
in both branches, which is in general different from the 
mass parameter in in the action. This provides a proof of 
our previous claim that the theory is not overconstrained. 
Thus, not only the background equation of motion but 
also the quadratic action for tensor perturbations are ex¬ 
actly the same as those in the dRGT theory Hi- On 
the other hand, for scalar and vector perturbations, there 
is no additional modes stemming from the gravity sector. 
This is consistent with what we have found in the previ¬ 
ous section, namely the fact that the number of physical 
degrees of freedom in the gravity sector is only two. (We 
shall discuss each branch later in this section.) Since the 
only two physical modes from the gravity sector coincide 
with the tensor modes, the cosmological background is 
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stable, provided that [i 2 > 0 and that the matter sector 
is stable (except for the standard Jeans instability that 
drives the structure formation). In particular, the the¬ 
ory automatically avoids the nonlinear ghost instability 
found in [hj and the classical Higuchi ghost @. This fea¬ 
ture of the minimal theory is in a sharp contrast to the 
dRGT theory. 

In the minimal theory, the constraints Co and Ci play 
key roles in eliminating the unwanted helicity-0 and 
helicity-1 gravitational modes. Actually, we have chosen 
Cq and Ci so that they always contain the two constraints 
C T (t = 1,2) in the precursor theory as well as two 
additional constraints. Moreover, as stressed just after 
m, Co ~ 0 and the Bianchi identity result in the same 
equation for the homogeneous and isotropic background. 
These rather non-trivial properties uniquely characterize 
Co and Ci. 

Phenomenology in the self-accelerating branch of the 
minimal theory is almost the same as the standard 
ACDM cosmology. For the background evolution, since 
X is set to be a constant by ED, the graviton mass term 
in m behaves as an effective cosmological constant, that 
can drive the acceleration of the cosmic expansion even 
without the genuine cosmological constant. Scalar and 
vector perturbations behave in exactly the same way as 
in the standard ACDM cosmology. Only the tensor per¬ 
turbations are modified by the graviton mass term, as 
described above. Therefore, the self-accelerating branch 
of the minimal theory leads to absolutely stable and phe¬ 
nomenologically viable cosmology. 

The normal branch of the minimal theory is also inter¬ 
esting. This branch is defined by X + NHX — MH = 0, 
and thus in this branch X is not a constant in general, 
making the contribution of the graviton mass term to 
the Friedmann equation (1121) dynamical. It is possible to 
find regimes of parameters in which /x 2 is positive and 
the effective equation-of-state parameter of the graviton 
mass term is either larger or smaller than — 1. This is 
quite remarkable both theoretically and observationally. 
It is commonly believed that the effective equation-of- 
state different from —1 would imply the existence of extra 
helicity-0 (and possibly helicity-1) degree(s) of freedom 
in either the dark sector or the gravity sector. Indeed, 
this is one of the main motivations for numerous dark 
energy surveys in the world. On the contrary, in the 
minimal theory of massive gravity there is no extra de¬ 
gree of freedom, and yet the cosmology with the effective 
equation-of-state different from —1 is possible and stable. 

In the absence of extra gravitational degrees of free¬ 
dom, in either branch of solutions, the minimal theory is 
not constrained by fifth force experiments and thus the 
bound on the graviton mass is relatively weak. We do 
not even need screen mechanisms such as Vainshtein’s 
one 0- The strongest bound, 

(15) 


comes from modification of the emission rate of the grav¬ 
itational waves from binary pulsars, which is of order 
O^l/u 2 ) 0. Here, /x s is the mass of gravitational 
waves at the source position and u> is the characteristic 
frequency of the system. 

The theory opens up new possibilities for gravitational 
phenomenology. One such example is possible appear¬ 
ance of a sharp peak in the stochastic gravitational wave 
spectrum 0. 

Conclusion. — We propose here a minimal theory of 
massive gravity, where only two physical modes in the 
gravity sector propagate, the gravitational waves, which 
become massive. This theory, contrary to the dRGT the¬ 
ory, allows stable homogeneous and isotropic (FLRW) 
cosmologies to exist. In particular, the self-accelerating 
cosmological solution, that was originally found in the 
dRGT theory, is now rendered stable in this theory at 
the fully nonlinear level. Furthermore, we expect that 
its phenomenology will be closer to General Relativity, 
as only the tensor modes are dynamical, as in General 
Relativity. The phenomenology of gravitational waves 
differs from the one in General Relativity, as they pos¬ 
sess a non-zero mass, giving them a different propagation 
dynamics (e.g. the propagation speed for modes with low 
frequencies). The normal branch cosmological solution 
is also rendered stable in the minimal theory and allows 
the effective equation-of-state parameter to differ from 
(either larger or smaller than) —1. As far as the authors 
know, this is the very first example in which the effective 
equation-of-state parameter is made different from — 1 
without introducing any extra degrees of freedom. This 
potentially has rather strong impact on various dark en¬ 
ergy surveys in the world. 

It is possible to extend this theory to a bigravity theory 
with only 4 physical degrees of freedom, by promoting 
M, M k E L k , and E L k to be dynamical. Two of the four 
are massless graviton degrees and the remaining two are 
massive graviton degrees. It is also straightforward to 
generalize it to a multi-gravity setup. We hope to extend 
this method to other cases for which unwanted degrees 
of freedom can be consistently removed via well-imposed 
constraints. 

Even within the context of single graviton, there are 
many possibilities for extension of the minimal theory 
proposed in the present paper. For example, one might 
be tempted to combine the minimal theory with the idea 
of the generalized massive gravity 0- Namely, by in¬ 
troducing Stiickelberg fields and promoting the constant 
coefficients 00 , 1 , 2 , 3,4 to functions of them, the dynamics of 
FLRW background can be altered. It is also possible to 
promote cp. i 2.3 .4 to functions of other dynamical scalar 
fields as in [14 Il5| . We leave the important discussion 
of the nature of the Lorentz violations, and their con¬ 
sequences related to the cutoff scale of the theory as a 
project to be discussed elsewhere. Here, we simply state 
that classically there is no Lorentz violation in the mat- 


/x s < 10“ 5 Hz, 
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ter sector and that Lorentz violation via loops should be 
suppressed by the tiny factor m 2 /Mp, provided that all 
matter fields couple minimally to the physical metric g 
at the classical level. 
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